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Abstract 



In this paper we give explicit formulae in momentum and coordinate space for the three- 
nucleon potentials due to p and tt meson exchange, derived from off-mass-shell meson-nucleon 
scattering amplitudes which are constrained by the symmetries of QCD and by the experimental 
data. Those potentials have already been applied to nuclear matter calculations. Here we display 
additional terms which appear to be the most important for nuclear structure. The potentials 
are decomposed in a way that separates the contributions of different physical mechanisms 
involved in the meson-nucleon amplitudes. The same type of decomposition is presented for the 
TT — TT TM force: the A, the chiral symmetry breaking and the nucleon pair terms are isolated. 



*Permanent address: Centro de Fisica Nuclear, 1699 Lisboa Codex, Portugal. 



1 Introduction 



Three-nucleon potentials based on tt and p meson exchange have been derived from different under- 
lying approaches. The recent papers which include p exchange [||, |2|, ^, ^ tend to extend to the rho 
meson whichever approach had been already used for the two pion exchange potential considered 
earlier (because of its longer range). The early history of the vr — vr three-nucleon potential was 
categorized in this way and summarized in [^. Those potentials built upon the excitation of the 
nucleon into a A isobar were the first and contain the least physical input. These A mediated 
potentials have recently been extended to include p exchange in Refs. and [Q]. Soon after the 
first A mediated tt — vr potential was constructed, it was realized that chiral symmetry (breaking) 
must be included in the potential. That is, the theoretical vrN amplitudes must satisfy well-defined 
chiral (so called "soft pion" ) limits |Q . One way to accomplish this begins with the use of (effective) 
Lagrangians with pions, nucleons, deltas etc. which satisfy approximate chiral symmetry ||^, |8|, |^. 
This approach has recently been extended to include p mesons, so that the assumed Lagrangians 
must also obey gauge invariance ||2|. This approach has built in the correct symmetries but the 
couplings must be estimated and, furthermore, the results have not really been tested against the 



pion-nucleon data [10|. 

Alternatively, the approach used in the Tucson-Melbourne (TM) family of forces is based upon 
applying the Ward identities of current algebra to axial- vector nucleon scattering. The Ward iden- 
tities are saturated with nucleon and A(1230) poles. Then employing PCAC (partial conservation 
of the axial- vector current), one can derive expressions for the on-mass-shell pion-nucleon scatter- 
ing amplitudes which map out satisfactorily the empirical coefficients of the Hohler subthreshold 



crossing symmetric expansion based on dispersion relations |11]. The off-mass-shell extrapolation 
(needed for the exchange of virtual, spacelike pions in a nuclear force diagram) of the most im- 
portant amplitude can be written in a form which depends on measured on-shell amplitudes 
only. This rewriting of the amplitude to bury all reference to models of the A exploits a convenient 
correspondence between the structure of the terms corresponding to spontaneously broken chiral 



symmetry and the structure of the model A term. However, both the field theoretic [12| and dis- 
persion theoretic ||13[ A contributions to result in an equally excellent description of the data 
|14| , so the early emphasis on maximum model independence can perhaps be relaxed. This will be 
done in the present paper. 

Extending this current algebra-PCAC program to the Compton-like processes |jl^ + N ^ 
+ A^, where is a vector photon and is an axial- vector current, provides, through vector 
dominance, a solid basis for modeling a p — vr exchange three-nucleon force. Again the amplitude 
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is obtained via Ward identities but this time exploiting in addition the gauge condition on the 
off-shell electroproduction amplitude. Unfortunately, there is no empirical subthreshold expansion 
of the invariant amplitudes of pion photo- and electroproduction. Instead, the amplitudes are 
tested against the soft pion theorems and (on-shell) against multipoles at threshold or slightly 



above threshold. The amplitude which forms the basis of the TM p — tt force meets these tests |15|, 
including the challenge of the recently measured neutral pion photoproduction yi6[ . 

The so-called Tucson-Melbourne (TM) family of three-nucleon forces, in the particular case of 
p — TT and p — p exchange, has never been displayed in a form suitable for applications in general 
nuclear systems (in particular, we have in mind few-nucleon systems). In the work of Ref. Q those 
potentials were applied to nuclear matter where certain approximations could be carried out. As a 
result, Ref. Q failed to provide a simultaneously complete and convenient explicit form of the TM 
p — TT and p — p forces, ready for immediate use in any calculation, as we propose to present in this 
paper. 

Nevertheless, the nuclear matter calculation of ref |Q was indicative of the role of the p meson 
in screening the vr exchange effect. The same effect had already been seen in the nuclear matter 
calculation of Ref. [Q] where a three-nucleon force constructed only from A^A^ — AA transition 
potentials was used. 

Furthermore, calculations of the triton bound state with inclusion of the TM 27r exchange 
potential yield over-binding and strong dependence on the irNN form factor regulator mass |17, 18|. 
This reflects the presence in the vr— vr TM force of a term (usually denoted in the published literature 
as the "c" term ) in part due to chiral symmetry breaking and responsible for the rapid variation 
of the lower-order terms in the vr — A amplitude, as demanded by the soft pion theorems. It appears 
then natural to extend the calculations on the three-nucleon system to the inclusion of the TM 
p — vr and p — p forces hoping to cancel part of the vr — vr force effect, in a similar fashion to what 
happens in the two-nucleon interactions. Moreover, since the leading terms of the expansions of 
the the pN — > vrA and pN pN amplitudes do not have to build in a drastic variation in the 
low energy region, they do not have such short range singularities. In this spirit, the Sendai group 
included a p — it force in a triton calculation |19| which, however, due to lack of a consistent body 
of information in the literature, mixes several terms and parameter prescriptions of more than one 
origin. 

Another motivation to have p — it and p — p forces derived from the same principles used in the 
construction of the TM vr — vr force is to have available a more complete two-meson-exchange three- 
nucleon interaction, which is needed to draw conclusions about the role played by the delta-isobar 
in the force. Only after introducing the consistent full family of forces in the hamiltonian, does 
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it become legitimate to make comparisons with the Hannover (two-body A^A^ NA transition 
potential) approach and judge how much effect comes from true A-isobar propagation. A first 
study of this question, restricted however to the use of the TM tt — vr force, appeared most recently 



in the literature |20|. Another one is in preparation where the p — n and p — p TM force input is 
needed |2^, |2^. In the context of these recent developments in few-body calculations, the present 
paper supplies the necessary information to address the one-decade unsolved problem of comparing 
and constrasting the Hannover with the Tucson-Melbourne approaches. 

Finally, and aiming beyond the triton bound state problem, the stage of nowadays powerful 
calculations allows the information of this paper to be applied to 1) the present concentrated effort 
on experimental searches for direct evidence of three-nucleon forces in the three-nucleon continuum 
guided by the three nucleon-scattering calculations of the Bochum group and 2) the a particle 



wavefunction calculation ||2J], under the strong indications from existing calculations [^, 2(:, 27, ^ 
that the three-nucleon force effects, when fully considered, will be relevant in bringing the binding 
energy closer to its experimental value. 



2 The TT — TT Potentials 

In this section we carry out the program of splitting the vr — vr force into the contributions relative 
to the different physical processes underlying the vrN amplitude. Historically the current algebra 
amplitudes were given in a form which emphasized their model independent character. That was 
achieved by not disentangling explicitely the A contribution from the chiral symmetry breaking 
terms. Such an attitude was justified by allowing a more direct relation to empirical quantities, 
where that separation could not be seen. The process of undoing this way of displaying the force 
is implicit in the set of initial papers |5|, |2^. Here we merely summarize the procedure by showing 
the main steps. 

The TT — vr force was derived from a non-relativistic reduction of the Feynman diagram of Fig. 
1. We begin by defining the T-matrix for the three-nucleon process depicted in Fig. 1 in terms of 
the S-matrix: 

where Nf = Y\j(l/2Ejy^'^ {2m)^f ^'^ (the index j labelling all the particles in the final state and Uf 
being the number of fermions) and Ni (defined similarly) are the normalization factors necessary to 
make Tfi the covariant T-matrix. Then the matrix elements of the momentum space three-nucleon 
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potential W are given by the non-relativistic reduction of the three body 5 — 1: 

< P'lP'2P'3\S - l|PlP2P3 >NR= -i(2vr)(5(p'io -P20 -PSO -PlO -P20 -P30) < P' iP' 2P' 3\W\piP2P3 > , 

(2) 

so that 

< p[,P2,P3\W\pi,P2,P3 >= (27r)^5^(p; + p'2 +P3- Pi - P2 - P3) < P'iP'2P'3\Tnr\piP2P3 > ■ (3) 

We note that this choice for the normahzation corresponds to the use of Dirac spinors normahzed 
as ■uii = 1 and to momentum eigenstates such that < p\p' >= {2Tr)^6{p — p'). Most momentum 
space three-body codes normahze the momentum space eigenstates to the Dirac delta function 
without the factor of (27r)'^; this point is discussed rather extensively in Ref. ||2^. Our choice 
of convention does not affect the coordinate space formulae as it is (correctly) absorbed into the 
Fourier transform, such that < r\r' >= 6{f— r'). The full three- nucleon potential W is the sum of 
three cyclic terms Wi, W2, and W3 where W3 has particle 3 in the middle (active nucleon) as in 
Figure 1. The other two terms are obtained by cyclic permutation of the indices of the incoming 
nucleons. 

Now we turn our attention to the intermediate state depicted as a blob in Fig.l. It is described 
by the (two-body) vrN scattering amplitude ir^iq) + N{p) — > iT^{q') + N{p') written in terms of 
invariant amplitudes F and B which have the the most general isospin decomposition: 

= -u{p',s'){5'^{F+ - ^[l4'])+ie^^>^{F- - ^[^, /])}u(p, .) (4) 

4m 4m 

To proceed, one removes the nucleon pole from the four off-shell crossing symmetric invariants 
F^{y,t,q^ ^q'"^) and B^{i',t,q'^,q''^). Given the four background vr N amplitudes (obtained after 
the nucleon pole is removed) F^{iy,t,q'^,q''^) and B^{v,t,q'^ ,q''^) {v = {q + q') ■ {p + p') / {-im), 
t = {q — q')'^), one performs an expansion in the pion four momenta q,q'. In this expansion, after 
terms of the order of (/x/m)^ (/x and m stand for the pion and the nucleon mass respectively) or 
higher are dropped, only the (t-channel) isospin even non-spin flip F^ and the isospin odd spin-flip 
B^ amplitudes survive |2^. We will concentrate then on those two. The non-spin flip even 
current algebra amplitude is: 

F+iu,t,q'',q'^)=fiu,t,q'',q''')^ + C+{,^,t,q\q'^) (5) 

where a is the pion-nucleon a term and /t^ ~ 0.93 MeV. The double divergence q' • • q/ of the 
background axial vector amplitude denoted by contains the higher order A isobar contribution. 
This amplitude is given by ^ 
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,2, 2g*^ v^a{q' ■ q) 



9m'^ ^ 

- ^^^j^p^^[{M-m){M + mf + {q^ + q'^){2M + m) 

- q' ■ q{2M - m + {q^ + q'^){M + m)"^)] , (6) 

where g* f« 1.82/2"^ is the vrAN couphng, M = 8.825/i the A mass, m = 6.726/U the nucleoli mass, 
UA = {M-m'^-q'-q)/2m , a{q'-q) = [El-q' ■q)[{M+mfq' ■q] + {li'^-q' ■q)[{M+'m)E^-\{jx'^ + q' -q)] 
with i?7r = (-^^ ~ T^"^ + /U^)/2M the center of mass energy of the pion at the A resonance. Here we 
have taken pL^ = l/2{q^ + q'"^) which reduces to when both pions are on mass shelh 
In general, must have the simple form |13, ^ 

C+{u, t, q\ q'^) = + C2q ■ q' + 0((?^) . (7) 

On the other hand, the assumed form of the function /, 

t, q\ q'^) = (1 - - 1) + - 1) (8) 

(adapted |l4| for vrN scattering from the SU (3) generalization of the Weinberg low energy 
expansion for vrvr scattering [^]) is such that does satisfy the soft pion theorems 

F+(0,0,0,0) = , (9) 

F+(0, ;u2, 0, 1?) = F+(0, /i^ /i^ 0) = , (10) 

where C"*" vanishes identically, and (with the aid of Eq. (|7|)) the constraint at the (on-shell and 
measurable) Cheng-Dashen point: 

F+{f),2^l\^^^^^^) = ^ + 0{q^). (11) 

In contrast to meson-meson scattering, 13 in (^) is not determined by soft-meson theorems (be- 
cause the nucleon four momenta cannot be taken soft) and is to be extracted from experiment. The 
most recent data analysis from the Karslruhe group [|ll[ makes this extraction slightly dependent 
on t: 13 varies from 0.46 (for t = 0) to 0.52 (for t = fi^). This situation is different from the situation 



of older data 13, ^l| which was consistent with an almost constant /? value of 0.4 in the same 
range of t. 

Neglecting the u'^ and qo terms in (^ because they are of the order of (/x/m)^ or higher, the 
amplitude can be expanded in the three- vector pion momenta q and q' as follows: 

F^{0,t,q',q'') = -^ + (^^ - C2)q- q' - -^{q' + /) (12) 
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The last equation explicitly exhibits the separation between the (higher order in g^) A contri- 
bution — contained in the C2 term alone — and the remaining chiral symmetry breaking terms. In 
ref. |5| and subsequent discussions of the TM vr — vr force, the C2 and P constants in the coefficient 
of the q ■ q' term were eliminated in favor of the on-shell (measurable) quantity F~^(0, jj? , fj? , fj?) 

,2 



F-^(o,^^^^^^) = (l-/3)^ + ^ (13) 

Jit ^ 



According to Eq. ([T^), the A term alone does not bring a structure in momentum space different 
from the chiral symmetry breaking term in q-q'. To isolate the A contribution one simply evaluates 
the constant C2, which, from (^) (and taking both pions on-mass-shell), can be done by evaluating 
the derivative of (^) with respect to q ■ q', at t = {q-q' = 0) and u = 0. The result is 
C2 = — 1.54^~^. This coefficient coincides with twice evaluated at t = fi'^ {q - q' = /"^/2) and 
u = and compares well with the value — 1.58//~^, obtained for the total contribution from at 
t = {q • q' = l^"^) and v = {). The two last comparisons give an idea of the negligible importance 
of the O(g^) terms in Eq. (|^. 

From the vrN amplitude in conjunction with the vrNN vertices and pion propagators, one con- 
structs the three body T/v/j, which according to Eq. (^ defines the three-body force represented 
in Fig. 1 1^, The constant C2 contributes then to the overall coefficient "b" that has been used 
in nuclear calculations (6 = 60- + 6a; = C2) 

fc= +C2 = --^(^ -F+(0,/x^/.^;.2)) (14) 



The value b = — obtained from early experimental results |31] F^{0, fj? , /j,'^ , fj?) = —0.16/^^^ 
and (j/f^ = 1.13^^^ implies a chiral symmetry breaking coefficient b^j = b — C2 = — 1.04^~^. If the 
updated set of Karlsruhe data is taken instead, ^"^(0, /x2,/x2,/i2) = -0.28/u-^ and a/f^ = 1.03^"\ 
then b = —2.62^"^ and consequently ba- = —1.08/^"'^. The comparison between the numbers 6a 
and 6(j exhibits that the A terms do not dominate the amplitude and any description which does 



not obey chiral symmetry breaking |32| is insufficient. 

The isolation of the A from the broken chiral symmetry terms being complete, we have to turn 
to the backward-propagating nucleon Born term that is added to the background amplitude, to 
generate the three-body force. That term is calculated by subtracting the forward time-ordered 
propagating nucleon term from the covariant nucleon poles. Such a subtraction was done in refs. 



2£]. Dropping again terms that are of the order of (fi/m) , the backward-propagating nucleon 



term ( "Z-graph" ) is 



^2 
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where g is the pseudo-scalar ttNN couphng constant. Since F'^ shows the same momentum depen- 
dence as the third term in Eq. (|12|), the constant = —g^ /Am? = — 0.15^~^ contributes, together 



with the coefficient of the last term of Eq. (p^), to the overall constant "c" (c = + Cz) that 
multiplies the q + q' term, in the definition of the TM vr — vr three-body force: 

+F'^mn{^)^ (16) 



We note that the term proportional to F^j^^{0) did not appear before in Eq. (12). This term 
nevertheless is inserted in c because both the backward propagating part of the nucleon pole F^ 
and the A couple with the pion with a (assumed the same) form factor F^^^jsi^q"^) which is defined 
g{q^) = gF-j^MNio^)- The chiral breaking a term is a c- number [14| and has no intrinsic q^ 
dependence (although it is multiplied by f{v,t,q^,q''^)). It is convenient, if not necessary, however, 
since part of the amplitude is due to and to multiply the final amplitude by form factors, 
dependent upon q^ and q'^. Consequently, the constant term {a / f-,^"^ , labeled "a" in the literature) 
attains a spurious momentum dependence from the form factors. The term proportional to i^^^^(O) 
in Eq. (^) is inserted to correct for this spurious momentum dependence to the orders in q^ and 
q'^ kept in the amplitude. 

The value c = Ca + Cz = 1.0/i^'^ that has been used before changes slightly to c = 0.91/i~^, 
with the new determination of cr/ f"^- From these numbers one concludes that in the non-spin- 
fiip amplitude amplitude the chiral symmetry breaking term dominates the relativistic backward- 
propagating nucleon term. 

The spin-flip current algebra amplitude B~ is simpler to decompose into the physical contribu- 
tions. This amplitude is |13, ^] 

B-{v,t,q\^^) = {FY{t) + F^{t))^ - -l^+D-{v,t,q\q'^) , (17) 

where the function contains the higher order A contribution and is given in references ^[l^, |29| . 
After the non-relativistic reduction of the vrN T matrix is done, in the process of doing the fi/m 
expansion to derive the three-body-force, we note that B~ is multiplied by which is already 

of second order in fx/m, so we need only keep the first term in the expansion of (|l7[). Therefore, 
the electromagnetic isovector form factors are approximated by a constant value 

Fi^ + Fj^ = l. + 3.70 + O(^) (18) 

and only the constant term {Dq = 4.87^^^) in the expansion of D~ is kept. Taking into account 
the multiplicative factor —l/2m from Eq. (0), the A contributes with a coefficient d/\ = Do/2m = 
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0.36^ 



^. The remaining contribution, due to the two first terms of Eq. (JITD, is d, 



0.2 V _ 



The sum = dca + c^a = 0.60^~^ defines the ^4 parameter that one encounters in the published 
equations for the force. 

Finally the contribution from the backward propagating nucleon pole ("Z-graph") is given 



Traditionally then these disparate contributions have been added into a total "d" coefficient, d = 
di+dz, where d^ has already been defined and dz = —B~^ /2m singles out the pair term contribution. 
The constant dz has the same value as the pair term coefficient dz = Cz = — 0.15/x~^ of the 
amplitude. The total "d" term coefficient becomes d = dca + d/^ + dz = — 0.75;U~'^. Again the 
numbers indicate that the pairing term is small relative to the background amplitude and although 
the A is important, it is only half of the "d" coefficient. 

Now that we have undone the usual representation of the Tucson-Melbourne vr — vr three-body 
force, let us close by recalling the form of the amplitude vrA^ — > vrA^, 



which best expresses the degree of model dependence of the amplitudes. Here Tz stands for the 
"Z-graph" contribution, Tz = Tb — TppB, where Tb is the covariant Born term and TppB is the 
forward propagating (positive energy) nucleon term. The model independent AT is added to the 
full Born term Tb so that to leading order Tb + AT satisfies the low energy theorems (Ward identity 
constraints) exemplified by (|9|) and (|lO|), and q' ■ M ■ q is a. background term for which (isobar) 
models are necessary. The grouping of the isobar contributions into AT + q' ■ M ■ q in Eq.( pO|) 
enforces the largest degree of model independence of the off-shell amplitude. This is because to 
lowest order (all) isobar contributions are included in AT, and q' ■ M ■ q is constructed to give 
contributions only in higher order terms. These general statements apply equally to the p — n and 
the p — p amplitudes discussed later on. 

What is unique about the ttN T"*" amplitude is that the background term q' ■ M -q oi Eq. (Q) has 
just the right structure to combine nicely with the momentum dependence (§) of AT so that the 
entire off-shell amplitude can be obtained directly from data without a specific reference to the A 
contributions to either AT or to q' ■ M -q. This point, while giving us greater confidence in the three- 
body force derivation, has lead to considerable frustration among those interested in three-nucleon 
forces. Hence the deconstruction above which is summarized in the following table. The table 
organizes vertically the coefficients according to the type of momentum and spin-isospin structure 



to zeroth order Q by 




(19) 



T = Tz + AT + q' ■ M -q 



(20) 



9 



which ultimately they multiply. Within each row, the decomposition acording to ([20|), is provided. 
The model independent terms labeled AT in the table arise from current algebra but are of two 
different types. In the first three rows (labeled a,b,c) they contain the pion-nucleon sigma term, an 
axial current-axial charge commutator which is a measure of chiral symmetry breaking. The last 
row (labeled d and originated by B~) contains a AT which is also a current-current commutator. 
However in this case the algebra is closed by the electromagnetic form factor of the nucleon and 
has nothing to do with chiral symmetry breaking. As we shall see in subsequent sections, the 
chiral symmetry breaking contributions to the AT of the pN ttN amplitute are quite small 
and, of course, absent altogether in the pN — > pN amplitude obtained via vector dominance from 

The specific momentum and spin-isospin structure of the final three-nucleon force is provided 
in momentum space in Ref. ||2^ and coordinate space in Ref. |^3|, and will not be repeated here. 
The strength constants of the 2tt exchange three-body force are given in units of the charged pion 
mass (139.6 MeV). The potentials, however, use an isospin formalism instead of charge states so 
it would seem natural to employ the SU(2) average pion mass (2m^+ + m^o)/3 = 138 MeV in the 
propagators and form factors. 
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3 The p — Ti Potentials 

The working definition of the p — ir three-nucleon potential is, as in the vr — vr case, Eq. (|3|) written 
in section 2. The three-body transition matrix T^^ pictured in figure 2 contains now two different 
meson propagators together with the vrNN and the pNN vertices from the external nucleonic lines of 
the diagram, in conjunction with the pN ttN intermediate amplitude. The amplitude is modeled 
through vector dominance from the amplitude — > vrA^ which is obtained by the method of Ward 
identities |^]. Those identities are (single and double) divergence conditions that determine the 
terms of the amplitude that are constant or linear in the exchanged momentum. Of the several 
spin-isospin component of the p — tt three nucleon potential only a few appear to be numerically 
significant. They correspond to the leading terms of the low energy expansion of pN — > vrA'^ 
truncated to terms of second order, as was the vr — vr force, in powers of the exchanged three- 
momenta between the nucleons. The (t-channel) isotopic spin decomposition of the amplitude is 
the same as that of the vrN amplitude of Eq. @ because the p is an isotriplet as well as the 
pion. Again, the expansion can be rearranged into the form of Eq. ( [20| ) for an optimal display of 
its degree of model independence. The isospin even model-independent terms + AT'^ nearly 
cancel to leading order because AT"*" goes to the soft pion — > 0) p analog Fubini- Furlan-Rosetti 



limit of pion photoproduction [15, 16|. The total model independent term in the A; ^ , q^O 
limit takes on the form 1 — mgA{0) / (fng) which is of the order of few percent (Goldberger-Treiman 
discrepancy[^). This cancellation is lessened away from the expansion point at zero, but then the 
model-dependent (and A dominated) background term q' ■ M ■ q appears to totally dominate the 
pN ttN amplitude. Therefore, in the isospin even amplitude, only the A-contribution, obtained 
in a dispersive sense, is developed into a p — tt potential in this paper. The situation is just the 
opposite for the model independent terms + AT^ of the isospin odd amplitude. Here the low 
energy theorem of pion-photoproduction is due to Kroll and Ruderman and is simply the "Z-graph" 

in pseudoscalar coupling This leading order (constant) term from the the p-analog Kroll- 
Ruderman theorem as g — > , /c — > should then form the basis of an important part of the p — tt 
potential. It must be supplemented, however, by the largest /) — vr potential terms arising from the 
A-dominated isospin odd background. Potential terms other than the above from AT^, including 
the chiral symmetry breaking terms, appear to be much smaller than these just mentioned and will 
be neglected in the following (for details see Ref. [^). 

We now pick out the terms with one p and one vr from Figure 2a which are believed to dom- 
inate the expansion and write them in detail. Those terms correspond to the "Z-graph" and the 
seemingly largest terms from the A isobar. The contribution from Figure 2b follows from 2a with 
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the substitutions 1^2 and k ^ —q where k is the three-momentum of the rho and q the three- 
momentum of the pion. (In what fohows the notation is shghtly different from that of ref. Q but 
closer to existing codes and practices). 

3.1 "Kroll-Ruderman" p — tt potential 

The analogue of the Kroll-Ruderman term of pion-photoproduction forms the basis of a specific 
local contribution to the overall p — n exchange three-body potential that presents a first order 
dependence on the exchanged mesons momentum k and q. It is consequently a contribution that 
is expected to compete strongly with the delta-isobar terms that are quadratic in the momentum. 

This contribution comes from the t-channel isospin odd term in the pN — > vrA^ amplitude. When 
the pseudoscalar coupling for the irNN vertex is assumed, this term is the backward-propagating 
nucleon Born term. In a form that evidences sequentially from left to right the role of the four 
vertices involved in the diagram, we can write for its momentum space representation 



_ 3 5^{p[ + p'2 +P3- Pi -P2- fs) 

{q^ ^ p^){k^ ^ mV' 



<p'i,P2,P3lWp"^(3^)|pi'P2,P3 > = +(2^) 



pi 

i f{ ■ f2 X 1% idi X k ■ (73(T2 • q 

F^,{k') F2,{k^))F^,{e)g'F^^^{q') (21) 

From this form it is evident that the coupling of the p meson to the inner "blob" is proportional 
to Fip{k'^) , in accordance with the low energy theorem of Kroll and Ruderman. Also, to obtain 
Eq. (^) we kept the distinction between the Dirac and the Pauli couplings of the p meson to the 
outer nucleon, which appears as the sum Fip{k'^) + Kp F2p{k'^) ■ Subsequently, we retain in all 
equations one form-factor Fip{k'^) for the Dirac coupling of the rho to the nucleon and another 
form-factor F2p{k'^) for the Pauli coupling. The latter is multiplied by the on-mass-shell Kp which 



is taken to have the Hohler and Peitarinen value 6.6|35]. Separate form factors are not usual in 
nuclear physics, where it is common to take the same form for both the charge and the magnetic 
momentum form factors. Our choice has the advantage of being general and fiexible to any type of 
(independent) behavior of the two form factors that experiment may eventually reveal. The price 
we have to pay is the introduction of a new parameter relative to F2p{k'^) . 

The coupling of the pion to the nucleon is described by a form- factor Fj^n . Other form factors 
are introduced as needed. All form factors are normalized to unity at the meson on-mass-shell 
momentum. As for the structure, (T2 • g is the non-relativistic coupling of the vr to one of the outer 
nucleon lines, while di x k comes from the p coupling to the other external nucleon. 
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In coordinate space the "Kroll-Ruderman" potential becomes 



g^gi rl • X f2,{(fl X Vl • (T3)C72 • V2 



(47r)2l6m3^ ""p 

{{ZoDpixi^) + KpZpDp{xi3))Zf,{x23) 



(22) 



where Xij = x 



Xj and Vj = V^3 



d/dxl^ = j^^d/dxis, {j = 1,2,3 specifies a given cartesian 
component and i = 1,2 a given miclcon). 

The gcncrahzed form factor H{q^) in each dimensionless coordinate space function Zp{xij) is 
typically a product of monopolc form factors: one for the vertex to the outer nuclcon and a second is 
included in the coupling of the meson to the "blob" of figure 2. The coordinate space "propagator" 
has the generic form 



47r r d^q H„/3{q^) 
IX J (27r)3 g2+^2 



(23) 



where 



Haf}{q' 



(24) 



For the exchanged pion Ha/3{q^) = F^j^^{q^) so that = Kp and the a,/3 indices are suppressed. 
The p propagator with subscript "DDp" also indicates Aq, = A^ with values from the Dirac coupling 
( Fipik"^) )2 and the transformation /x — rup but the subscript "PDp" indicates Aq, 7^ A^ with values 
arising from the product of F2p{k'^) and Fip{k^) and /x — rUp. 

With these form factors and carrying out the indicated derivatives, the "Kroll-Ruderman" 
potential becomes 



< r[,r'2,r'3\W-{3A)\f{,r2,r3 > = 



^"^g'gpH^^ - xi)dHxl - rr-2)<5=^(4 - X3) 



(47r)2l6m3^ 

n • T2 X 7^(crl X Xi3 • a3)a2 ■ X23 
((^lV(^13) + «p^i%(^13))Xr(x23) 



where Xi{x) = Z'{x). Explicitly for A^ = A/3 = A 



and for Aq, 7^ A^ 



xf (x) = -11 



Ginx) - { ^ 



G{nx) 



^C(Ax)-i(^-l)e-(-) 



2 / A? - ■•2- 



(25) 



(26) 



(27) 
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where 

G{x) = ^ (l + ly (28) 

Note that the generic X\{x) = Z'{x) has dimensions of mass and the spin operators are coupled 
with unit radial vectors. Then the overall constant has units of MeV and the magnitude is 28.72 
MeV. It should be clear that for rho exchange the generic X^{x) is written as /i — > nip. This finishes 
the explicit display of the "Kroll-Ruderman" term. It is expected to be the most important of the 
model-independent parts of the rho-pi potential. 

3.2 piT Potential with A intermediate state 
3.2.1 Isospin even amplitude 

Of the eight spin functions the potential decomposes into, the terms proportional to 1 + Hp are 
believed to be the largest because Kp has the value 6.6 on the rho-mass-shell. Those terms which 
are (t-channel) isospin even are generated from the spin functions ffi • q(T2 ■ Q and ffi • k(f2 ■ q. They 
take the momentum space form (see Eq. 2.19a of ref. 

<Pl,P2^P3\^pnAi^MPl,P2,P3> = -(2vr) — 



fi- f2{k ■ kai ■ qa2 ■ q - k ■ qai ■ k(j2 ■ q) 
■^9p{ Fip{k^) + Kp F2p{k^) ) Gl,^Fpr,A{k') 

m5M-rn^ ^^^^^^ 
M M — m 

where M is the mass of the Delta and m is the nucleon mass. The two terms can be combined 
as ((fi X k) ■ {q X k)(f2 ■ q which manifests the pNN coupling to "outer" nucleon 1 and the vrNN 



coupling to "outer" nucleon 2. On the other hand, we note that the second term of Eq. 29 has the 
same structure (ffi • k(j2 ■ qk ■ q ) as the "b" term of the 27r exchange potential. 
In coordinate space the potential becomes 



< r'l, r^, r^|W+^^(3^)|rl, ra, rg > = 5\x'^ - xi)5\x'^ - X2)5\x'^ - £3) 



n ■ 7^[Vi • Vidi ■ V2cr2 • V2 — Vi • VacTi • Vicfa • V2] 

{{ZdGp{xi3) + KpZpGp{xi3))ZN^^{x23) (30) 
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where subscripts "DGp" refer to the generahzed form factor arising from the product of Fip{k'^) and 
the ^AA^ form factor FpNA{k^) , ''PGp'' to the analogous product F2p{k'^) FpNAik"^) , and "iVA^" 
to F^NNiq^) F^NAiq^)- 

The second derivatives of Eq. ( pO] ) are most conveniently carried out with the aid of the following 
basic identity 

~ - ....... (3^) 



didjZ{x) = -[5ijY2{x) + {3xiXj - 6ij)X2{x)] 



{xi = Xi \xi\) where X2{x) = Z"{x) — l/xZ'{x), and 12(2^) = Z"{x) + 2/xZ'{x). Explicitly they are 
X^{x) 



Y,^{x) 



1 A 

F{px) - -^F{Ax) - i-Ax(^ - l)G{Ax) 



-fJ,X 



A3 



-Ax 



fix 



1 A , A2 2 ,e-^^ 

Ax(^ - 1 (1 - — )-; 

p^ Ax 2fl >2 A j^^J 



Because 



Fx = 1 + - + - 

X X X'^ 



(32) 



(33) 



for the seldom needed Aq, = A/3 = A we can make the identifications Xl^lx) = fi'^T{x) and 5^^(x) = 
fi'^Y{x) familiar from the Argonne-Urbana two nucleon force expressions ||3^. In the present case, 
however, Aq 7^ A^ for each propagator and we must employ the longer forms 



F(/ix) 3 
e-M^ A3 ^ A?, 



W 5 A 2 ^ 



, A3 A|-/i 



-A„x 



+ 



//3^A2 - A2^ A^x fi^^Al-Ay Ac,x 



/9 



(34) 



Before displaying the final potential in r-space let us introduce the "tilde" notation for the 
combination 1 + Kp which appears so often. For example, the p propagator in ( pO| ) which has two 
form factors in the pNN coupling and only one in the pAN, takes on the compact form 

ZNGpix) = ZoGpix) + KpZpGpix). (35) 

Since differentiation is a linear operator we can write, for example, 

X^{x) ^ X^j^cix) = XP^^ix) + KpXPp^ix). (36) 

to shorten the following equations. 
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Using Eq. (pi) the first term of Eq. (pOt) (from cfi • €[02 ■ qk ■ k ) becomes 



< r[,r'^, r-^l W+^a(3^)I^"I' ^2, r"5 >^ = 5^{x{ - xi)5^{x'^ - X2)5^{x'^ - x^) 

-f^mp 5M-m 

r-{-r-2[^^Y^{x,3)Y^^{x23) 



(37) 



The second term of Eq. ( pjO| ) has the same structure [dl ■ k(T2 ■ qk ■ q) as the "b" term of the 27r 
exchange potential and is 



6-^{x[ - xi)d'\x2 - X2)d'\x'^ - X3) 



5M -m 



Glipimg*) 



n ■ f2-[ ai ■ (f2Y^{xi3)Y^{x23) 

+ 5i2(xi3)X2^(xi3)^f (^23) + 5l2(£23)1^2''(^13)^2''(^23) 
+ (9(Ti • f 130-2 • £23^13 • ^23 - 5'i2(xi3) - S'l2(x23) 

- ai ■ a2)XP{xi3)X^{x23)] (38) 

The sum of these two terms agrees with the corresponding equation (2.23) of Identification 
of Eq. ( pSD with the coordinate space "b" term as given in ref. can be made by remembering 
that 5'i2(xi3) = 3di • xiad^ • X13 — di • (T^ and noting that [l2(x)—X2(2;)]/3 = Z'(3;)/3;. The overall 
coordinate space constant (including the powers of meson masses in the definitions of (x) and 
1^2'' (x)) for both terms ( pT] ) and ( pS] ) is 524.28 MeV. This number is huge because the pNA coupling 
constant G^j^is (rather artifically) defined on the p mass shell but is in reality obtained from the 



experimental G^^(A;2 = 0) = 14.7. That is Glj^{k'^ = 0) 



A/7 1 



GllpFpNA 



one finds Gl 



'Mp - G*M^ik^ = where A = 5.8/x and ~ ^.up. j^hc piu^ui^u ^j^^^^ 

is employed in the subthreshold region for space- like momentum transfer, where it is much smaller. 

3.2.2 Isospin odd amplitude 

Those largest terms which are (t-channel) isospin odd multiply the spin functions of Eqs. (2.13f ) and 
(2.23g) in Ref. |Q: m^/Cg = i^^i x ^ • 0^3(^2 • Q and m^lCj = idi x k • gifa • k(J2 • q. The particular 
combination containing the delta-nucleon mass difference takes the form m^lCj — q ■ km^K^ = 



where A 



= 0). Solving for G^^, 
5.6/i. The product 0%^ FpNA{k'^) 
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idi X k ■ qa^ ■ kd'2 ■ q — q ■ kidi x k ■ (73(72 ■ q- This can be rewritten as i(fi x k ■ [((T3 x q) x k](72 ■ q to 
show the spin structure connected with the several vertices most clearly, or finally as — ^((73 • q x k(7i ■ 
k(72 • q — (Ti X (T3 • qik)'^ (72 ■ q)- The third form is convenient for the coordinate space manipulations 
and has appeared before from non-relativistic transition potential derivations ||l], We will quote 
the momentum space potential with the last form of the spin functions 



{q"^ + ^i'^){k'^ + m?p) 
i[(73 ■ q X k(7i ■ k(72 ■ q- (7i x (73 ■ q{k)'^(72 ■ q\ 



1717 gF^NN{q^){m g*F^NAiq^) 39 

M M — m 

One can readily see that the first term has the same spin structure (73 ■ qx k(7i ■ k(72 ■ q as the 
■'d" term of the 2tt exchange potential. In coordinate space the isospin odd potential becomes 



[(73 • V2 X Vicxi • V1CJ2 • V2 — Cfl X (T3 • V2V1 • Vl(72 ■ V2] 
{{ZdGp{x13) + KpZpGp{xi3))ZN Ap{x23) (40) 

After the derivatives are carried out, the first term of Eq. (^0|), analogous to the "d" term of 
the 2tt exchange potential, takes the form: 

< ^^1, ^1, ^ll W-^A(3^)k"l' ^"^2, r1 >^ = (^^(4 - x-'i)d\sp2 - X2)-5^(4 - xli) 

firup M + m ^ 

f{-f2X f3[{(73 ■ X23 X Xi3)(cfi • Xi3){(f2 ■ X23) 

X^{xn)X^ix23) 

+ (01 • Xi3){(f2 ■ Xi3 X 0^3) X^(xi3)Xf (X23) 

2^23 

+ (f^2 • £23)(o^3 • ^23 X ol) X[(xi3)X^(x23) 

Xl3 

+((f2 • ai X al)— Xf (X13)— (X23)] (41) 

Xl3 X23 
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The second term of Eq. ( ^0[ ) does not have a counterpart in the 2tt exchange potential but 
instead takes the simple form: 



finip M + m ^ 

n • 7^ X f3[+{(f2 ■ X23)('73 " 0\ X X23)Y.{' {xi3)X^ {X23) 

+{a2 ■ a-3 X al)— y/(xi3)Xf (X23)] (42) 

3^23 

The overall numerical constant in both ( |4l| ) and is 109 MeV, about a factor of five smaller 
than that of the isospin even term ( for this reason this term was not included in the nuclear matter 
calculation of ref. in a trinucleon calculation however it may turn out to be meaningful). A 
model three-nucleon potential derived from the transition potential approach would have an exact 
ratio of 4:1 for the isospin even versus isospin odd contributions. 

One local term proportional to i[cJi x k-(j3(j2 ■ q\k-k remains from the expansion of the covariant 
expression for the amplitude with an intermediate A. The overall numerical constant is about a 



factor of ten smaller than that of (41) and ( |4^ ) so we will not consider it further. It requires a third 
derivative of Z or first derivative of Y2 in coordinate space anyway, which eventually suppresses it 
further. 

The total rho-pi contribution is obtained from VV""^(3A) by first adding the contribution from 
W^^{3B) = {W^'^{3A) : 1 ^ 2) and then taking cyclic permutations of three "active" nucleons . 



18 



4 The p — p Potentials 

The definition of Eq. ^ is once more the starting point. The three-body Tjy/j for the Feynman 
diagram of figure 3 involves the pN pN amphtude, together with two p meson propagators 
and pNN vertices. The calculation of the latter amplitude is obtained by generalizing the Thirring 
theorem for Compton scattering to the case of p mesons, through the current-field identity, under the 
assumption of vector dominance. The longitudinal character of the p has however to be considered 



as well as the fact that it carries isospin. To do that one proceeds as did Beg [|36| in his analysis 
of Compton scattering of isovector photons. The dominant term in the background (free of the 
nucleon pole term) Beg amplitude arises from the t channel p pole (3/9 coupling) and from a direct 
ppfifi contact (seagull) term. Those contributions constitute the AT part of Eq. (]20|). The A term, 
contributing to the quadratic terms in the momenta, cannot be fixed by any divergence condition 
and is introduced, explicitely separated from the remaining amplitude, in a model dependent way. 
The following subsections will turn to the detailed form of all these components of the pN amplitude 
of different (physical) origin. 

4.1 Model-independent potentials from low-energy theorems 

The apparent dominant terms in the model-independent parts of this potential are odd in the t- 
channel isospin. The term displayed below results from the coupling of transverse (T) p-mesons 
to the low-energy amplitude constructed by Beg and has been labeled the "Beg" potential in the 
previous nuclear matter calculations. Similarly to the Kroll-Ruderman term, the "Beg" term is 
linear in each of the meson momentum, and it turns out to be important when comparing to the 
A isobar contributions. In momentum space it takes the form 



< p'vP2,P3\Wtt\pi^P2,P3 > = +(27r) -jr- ri • ra X rs 

(k ^ -|- m-i)[K' ^ + vrtp) 



[(ai X k) ■ {a2 X k') X aa] F^p{0){l + 
^( F,,{e) + F2p{e) )( Fi,(fe'2) + F2p{k'^) 



(43) 



where F2v{^) = k-v = 3.7 is the isovector anomalous moment of the nucleon. While the term 
independent of F2V in -f'ip(0)(l + yr^) arises from the backward-propagating nucleon pole term 
(or pair term), the term in F2v{0) represents a direct ppNN contact (seagull) term in the underlying 
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pN pN amplitude. In principal the second term should be distinguished by a form factor of 
a different character {Fip{{k — k')'^) than the previous terms of the pir potential. In practice, 
this component of the pp potential has so far been calculated in the approximation of no form 
factors on the "active" nucleon and therefore includes only the pNN form factor at the external 
nucleons. This approximation is similar to the approximation made in the "d" term of the 2tt 
force. In that case, only the constant term of the current commutator expansion Fiv{t) + F2vit) ~ 
Fiv{0) + F2v{0) + 0{q^) « 1 + 3.70 was retained (see Eq. (18)). 

In coordinate space the "Beg" potential has a slightly different form from the previous ones 
displayed, because the propagator Z of Eq. has only a single monopole in the H function of Eq. 
(|23) and (as in the "Kroll-Ruderman" term) only the first derivative is needed. With the definitions 



A2 

G{px) ^G(Ax) 

p^ 



(44) 



where the added subscript 1 is meant to remind about the single formfactor, the coordinate space 
potential becomes 



or 



< ^1, rl, r|| W;|T,^|r1, ri, rg > = 6^{x[ - xi)5^{x'^ - X2)5^{x'.^ - x^) 

^^^^p^5pi^i,(0)(l + -p^) 

n • T2 X T3(crl X Vi) • {(T2 X V2) X 

Zlnip{xiz)Zlmp{x2z) (45) 



5^{x[ - xi)5'^{x'2 - X2)5\x'^ - xl) 

n • T2 X 75((fi X X13) • {(72 X X23) X (J^ 
Xfi(xi3)Xfi(x23) (46) 

4.2 pp Potential with A intermediate state 

There are two local potentials which correspond to the coupling of transverse p-mesons to the A 
pole. We introduce the notation 
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for the low-momentum transfer pAN vertex. The two local terms correspond to t-channel isospin 
even and odd amplitudes. We take them up in turn. 

4.2.1 Isospin even amplitude 

The form of the term generated by this amplitude agrees with the form arising from a transition 
potential construction and reads in momentum space 



< K , V'l , P3 1 l^^l ,P2,m> = - (27r) 



3 {p[ + p'2 + p's - Pl -P2- Ps) 

Ti ■ f2 {{d\ X k) xk) ■ {{(72 X k') X k') 

^^fu-m) ^ ^'^^^'^ + '''' ^"'^^'^ )FpNA{k') 

{ Fi,{k'^) + F2p{k'^) )F,MA{k'^) (48) 



which, after Fourier transform becomes 



< rl,rl,rl\wl'^l\f{,f*2, ra > = - x-{)6^{x^2 - x~'2)5\a/3 - xli) 

(47r)2 18(M - m) ■ 
[ {{ai X Vi) X Vi) . ((^2 X V2) X V2)] 

{ZdGp{x13) + KpZpGp{xi3)){ZDGp{x23) + KpZpGp{x23))-{^9) 

Carrying out the derivatives one obtains the final form of the matrix element in r-space 



< r[, r^, r-^|W^+i|rl, r2, ri > = 6^{x[ - xl)(5='(x^ - X2)(53(x^ - ^3) 

(47r)2 18(M-m) ^^"^^9 
([4 oi • a2Yi{xis)Y^'{x23) - 2 Su{xi3)X^{xi3)Yi{x23) 

-2 5'i2(x23)^2'(^13)-^2 (^23)] + ■ X23(^l ' ^130^2 • ^23 

- 5l2(xi3) - 5i2(f23) - Cn ■ a2]X^{xi3)X^{x23)]). (50) 
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4.2.2 Isospin odd amplitude 

There is a relative factor of —1/4 between isospin odd and isospin even in this case. This is precisely 
the relative factor obtained from a transition potential approach. The reason why the TM force 
ends up with the same factor lies in the fact that the k ■ C ■ k' contributions to the pN — pN 
amplitude were modeled by A-poles in the s- and u-channels. 



< P'l ,P'2,P3\ %pA \Pi ,P2,P3> = + (27r) 



3 5^ {p[ +P2+P3- Pi -P2- Ps) 



^ fi • tI X 7^ 0*3 • [ ((al X k) X k) X ((cfi x k') x k')] 

( + F^pik'^) )FpMA{k''') (51) 



In coordinate space the equation becomes: 



(47r)2 18(M-m)4 ' ^ ^ 

as ■ [ (((tI X Vi) X Vi) X {{(72 X V2) X V2)] 

{ZdGp{x13) + l^pZpGp{xiz)){ZDGp{X23) + l^pZpGp{X23)){^'^) 

We note that the expansion of the vector products in the previous equation leads to four terms, 
one of which is identical to the "d" term of the vr — tt force. 
After carrying out the derivatives we find finally: 



< r'^y^^ ?'3l>VipAkl' ^3 > = <5^K - Xi)S\x'2 - X2)S^{x'^ - xs) 



9p^M 



Tl • r2 X T3 



(47r)2 18(M - m) 4 
[(^2 • (Ts X a-{)Y2p{xi3)Yi{x23) 

-{as ■ (fi X X23){(T2 ■ A23^2'(3^13)^2 (3^23) 

-{ai ■ (72 X as)Y^{xis)^^ 

X23 
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-((73 • Xi3 X a2){(Tl ■ X13)X^{X13)Y.^{X23) 
-(a3-a-ixa2)^^y/(x23)] 

-[(0^3 • X23 X Xi3)(oi • Xi3){(f2 ■ £23)-^2 (^13)^2 (2^23) 

1 _ 

-(cr"l • ^13)(c^2 • ^13 X ffs) -^^2 (2^13)^1 (2^23) 

3^23 

-(0-2 • a:23)(o-i • 0-3 X 2:23) ^2(2^23) 

Xl3 

-(£72 • CJl X (Tsj (53) 

X\3 X23 

where the four last terms correspond to the "d" term analog piece of the matrix element. 

5 Recommended choice of parameters 

The set of numerical values for the coupling constants and the monopole form-factor regulator 
masses are compiled in Table II. The monopole formfactors are defined in Eqs. (23) and (24). 

This set differs from table A2 of ref. |Q only in the value for the vrNA coupling constant, which 
should be g* = 1.82/x~^, instead of the larger value of that table. Only this value is consistent with 
the vrNA coupling constant implicit in the 27r exchange three-body force (see Appendix A of Ref. 
1^]) and with the value obtained by the Karlsruhe group from ttN scattering analysis. They 
used a fixed-t dispersion relation for the invariant amplitudes in order to determine the parameters 
of an "effective pole" which simulates A-exchange at energies near and below threshold. The 
coupling constant so determined corresponds to the dimensionless coupling /* (/*^/47r = 0.26) 
used by Martzolff et al [p, and differs from the Hannover group choice f*"^ /Att = 0.35 for their 
transition potential The latter value is derived from the A- width, using the interaction term 
of the Lagrangian and "should not be used in application of pole term formulas near or below 



threshold" |38]. This latter value is 40% higher and therefore emphasizes the effects of the forward 
propagating isobar in the transition potential picture compared to the covariant isobar contribution 
to ttN scattering. 

The decision for the prescribed values of the remaining parameters was justified in ref. 
Briefly, we recall here that, 

• the coupling constants concerning the vr come from low energy ir — N scattering data. The 
value of (7 = 13.4 has been challenged in recent years by partial wave analyses of nucleon- 
nucleon scattering |3S]. Much discussion has ensued on this still unsettled point. 
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the cut-off value for the irNA form factor is taken to be the same as vrNN, which is suggested 
by data from neutrino-nucleon scattering (p + ^ A^"*" + ^~). 

for the p meson, the choices made on the parameters were determined through vector meson 
dominance models (which ties the needed hadronic pAN coupling to the 7AN vertex) that 
are believed to be valid for processes involving low four-momentum transfer. We follow the 
experimental analyses of the 7AN vertex and use a monopole form factor with cutoff mass 
ApATA = 5.8)U to approximate this vertex. As for the pNN vertex, it has two couplings: 
a direct Dirac coupling and the Pauli coupling to the anomalous magnetic moment of the 
nucleon. In the Tucson-Melbourne prescription, the ratio of the Pauli to Dirac coupling is 
allowed to vary with the four-momentum of the p. This choice appears to accommodate 
reasonably both the vector dominance value of = 0) = 3.7 and the on-mass-shell value 

of Kp{k'^ = vrip^) = 6.6 The two cutoff masses Ajjirac = 12/i and ApauH = 7.4// are 

then determined from nucleon electromagnetic form factor data. Somewhat later, a group led 



by Gari [40| made a detailed fit to the electromagnetic form factors of the nucleon with the 
ansatz suggested already in and came to a similar conclusion as . What was not noticed 
by the authors of |^ and presumably of |4C] was that the A;^ dependence of the Dirac and 



Pauli rho form factors away from the rho mass shell were also shown and even parameterized 
for spacelike momentum transfer by Hohler and Pietarinen in |35]. Indeed, the ratio of Pauli 



to Dirac coupling from the latter analysis does not vary much with k"^ and is 6.0 at /c^ = 



1 41]. It turns out that this form factor result for the pNN vertex is quite compatible with the 



data for the electromagnetic isovector nucleon form factors 

We next discuss the most varied aspect of the TM force: the fact that each user feels free to 
substitute his or her own choice for the recommended cut-off value of the monopole ttNN form 
factor. Its recommended value (A = 5.8p ~ 810 MeV) is inspired by the 5% Goldberger-Treiman 
discrepancy (which, measuring chiral symmetry breaking, is consistent with a 3% decrease in the 
ttNN form factor from = mj^ to = [^. The present understanding of the variation of 
the other quantities in the Goldberger-Treiman relation from the chiral limit to the chirally broken 
"real world" , imply an "error bar" of ±200 MeV for the A determination. This makes our choice 
in good agreement with the OPEP tail in many contemporary NN potentials: A ~ 980 MeV for 
Argonne V14 Q, A = 800 MeV in the present day evolution of the Bonn OBEP potentials Q, 
and A = 800 MeV for both vr and p in the Bochum potential The vrNN form factor A ~ 1300 
MeV of the full Bonn potential |45] has always been inconsistent with evidence other than partial 



wave analysis of the NN system this difficulty is expected to be overcome if additional diagrams 
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with correlated p — it exchange is included into the full model |4(:]. On the other hand, the cutoff 
mass Aq = 965 MeV in the exponential form factor of Nijmegen potentials |47] corresponds to 
A = \/2Ao ~ 1365 MeV for a monopole form factor |^], rather like that of the full Bonn potential. 

We should note that the current controversy over the value of the on-shell vrNN coupling 
g affects the cutoff mass determination from the Goldberger-Treiman discrepancy Aqt = 1 — 
''^9a{^) / {f-wd) ~ 0.05; a smaller coupling constant implies a smaller discrepancy and a larger cutoff 
mass. But other analyses which isolate one-pion exchange via a Regge analysis of the charge ex- 
change data da /dt\np minus da/dt\pp also find a 3% decrease in the ttNN form factor from q'^ = ttt-tt'^ 
to g'^ = 1^^, indicating a certain robustness in the recommended value of the cutoff mass. 

The recent consensus on a low mass cutoff for pion exchange highlights the point already 



emphasized by the Hokkaido group [25| and, in the modern context, by the Sao Paulo group |48|. 
The contact terms (those proportional to a coordinate-space (5-function and its derivatives) are 
spread out with increasing importance as A becomes smaller and the size of the nucleon grows. The 
dominant attitude so far was to consider that these contact terms, bringing the nucleon structure 
signature, should not be included in potential models. In contrast, the low cutoff mass of the vrNN 
and yoNN form factors of the Bochum NN potential, means that all mesons more massive than the 
p are subsumed into contact terms. This reformulation of the traditional OBEP can still yield a 
satisfactory description of the NN data, although a analysis is needed for a definite judgment. 
The nice feature of the Bochum NN potential is that the small two-pion exchange term is highly 
suppressed, cut down by four powers of a rapidly varying cutoff function and partially replaced by 
the contact terms. Weinberg obtains similar contact terms in his NN potential from a presumed 
four-fermion interaction Q. 

We finish this section by returning to the subject of the pNN formfactor. We have developed 
the formulas in momentum and coordinate space for two separate Dirac and Pauli formfactors of 
the monopole type. Clearly it is easy to specialize these general formulas to a single formfactor 
and then the choice of parameters is between the vector dominance value of Ky(k'^ = 0) = 3.7 
expected for a narrow resonance and the on-mass-shell value of Kp{k^ = nip'^) = 6.6 obtained from 
a dispersive analysis of vrN scattering. From our current understanding of the the pNN vertex, a 
better but more difficult to implement form factor might be those (Dirac and Pauli) of Eq. 4.3 in 
p5| which goes into a monopole for high k'^ but has a more complicated structure at small /c^ near 
0. We suggest that future users of three-body forces might concentrate on the pNN vertex rather 
than on the vrNN vertex, where the form factor of the Tucson-Melbourne three-nucleon potentials 
is not only determined by particle and nucleon-nucleon data but, in addition, is rather consistent 
with most of the realistic two- nucleon potentials. 
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6 Numerical Results 



It is well-known that calculations with many recent "realistic" two-nucleon potentials (Argonne V14, 
Nijmegen potential or the folded diagram version of the full Bonn potential) produce underbinding 
of the the trinucleon bound state [^]. However, the same interactions, when taken in conjunction 
with the TM two-pion exchange three-body force, provide too much binding. There are overlapping 
and retarded pion exchange graphs not discussed in this paper which yield non-local terms |^]. These 
have not been numerically evaluated yet, but are expected to be less important than the (already 
small) nucleon-antinucleon pair terms of the 27r-TM potential and henceforth expected not to alter 
the binding defect problem. 

To test the idea that the three-nucleon potentials of the two-meson exchange structure extended 
to include p-exchange act against the overbinding effect of the 27r-TM potential, we estimated 
perturbative contributions to the binding energy of the triton from the extended TM force described 
in detail in the previous sections. 

The calculation was done using a wavefunction in the coordinate-space obtained by solving the 



Faddeev equations with the Malfliet-Tjon I-III potential |50|. We present the results in Table III, 
where the individual contributions from the different physical processes considered in the force are 
singled out. The contributions from the vr — vr force are slightly different from preliminary results 
presented before |49] because in this paper we use the expansion coefficients of Table I. [^] We 
stress that the wavefunction used is "semi-realistic" at most (no NN p or d waves are included in 
the Malfliet-Tjon I-III potential). The past history of results obtained with the vr — tt force showed 
us already that any calculation with less than 18 three-body channels has to be considered with 
caution. The numbers given here were produced merely in the spirit of having a quick and hopefully 
qualitative information on the importance of the p exchange in the three-nucleon force. A much 
more sophisticated calculation, which is out of the scope of this paper whose aim is restricted to 



provide the force to be used, will be reported 22 1 soon. 

From table III we conclude first that the p exchange does indeed modify the effect of the 
TT — vr force and that the effect goes in the (expected) direction of less binding. Since it has been 
fashionable to take the derived strengths of a nuclear force seriously, but to consider the meson- 
nucleon form factor cutoffs as adjustable parameters, we show in the columns labeled "MartzolfF' 
a set of calculations which keep the strength constants of the Tucson-Melbourne three-body force 
but use the much heavier cutoff masses from Ref. [|l[] (we had to convert however their monopole 
parameterization of the product of two vertices to our dipole convention to find Ayr = 10.6// and 
Kp = 13.4/i for all couplings to N's or A's). We note then that although the contributions of the 
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TT — TT are, as is well known, strongly cutoff dependent, the inclusion of the vr-p force supresses, at 
least in this triton model, much of that cutoff dependence. Finally, the results obtained for the TM 
force show a satisfying pattern of decreasing effect with increased mass of the exchanged meson. 
The same pattern does not exist in the pointlike nucleons (large cutoff masses) calculation of the 
rightmost columns. In any event, the addition of three-body forces due to p exchange, especially 
the p — IT exchange force, does indeed counter the too strong attraction of the tt-tt and seems to 
stabilize the total effect under variations in form factors. We await future tests of these tentative 
conclusions with some eagerness. 
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Tz 


AT 


q'-M^-q 


Total 







+1.03 





+1.03 







-1.08 


-1.54 


-2.62 




-0.15 


+1.06 





+0.91 




-0.15 


-0.24 


-0.36 


-0.75 



Table 1: Expansion coefficients of the ttN amplitude used in the tt — tt force. Units of charged pion 
masses. Tz is the "Z-graph" contribution, AT is the model independent part needed to satisfy the 
low energy theorems, q' ■ M ■ q is the background term for which (isobar) models are necessary. 
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for the TT 

g* 1.8 /u-i 

^■kNN 5.8 
AttAta 5.8 /X 



for the p 



Qp 5.3 

Gmp* 191.1 

ApATAT (Dirac) 12 /x 

ApATAT (Pauh) 7.4 /x 

ApjvA 5.8 /X 



Table 2: Couplings labeled by a generic y and cut-off masses labeled by A. See Ref. 3 for complete 
definitions. 
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Exchanged mesons Tucson-Melbourne form factors "MartzolfF' form factors 



TT — 77 






AT 


q' -Ma-Q 


Total 




AT 


q' -MA-q 


Total 


r+ 

Total 


0.055 
-0.032 
0.023 


-0.409 
-0.052 
-0.461 


0.090 
-0.077 
0.013 


-0.264 
-0.161 
-0.425 


0.011 
-0.043 
-0.032 


-0.363 
-0.070 
-0.433 


-0.257 
-0.105 
-0.362 


-0.609 
-0.218 
-0.827 


p — 77 


T+ 
T- 
Total 


0.215 
0.215 




-0.022 
0.090 
0.068 


-0.022 
0.090 
0.283 


0.342 
0.342 




0.261 
0.340 
0.601 


0.261 
0.340 
0.943 


p-p 


T+ 
T- 
Total 


0.001 
0.001 


0.002 
0.002 


-0.005 
-0.007 
-0.012 


-0.005 
-0.004 
-0.009 


0.005 
0.005 


0.019 
0.019 


-0.101 
-0.247 
-0.348 


-0.101 
-0.223 
-0.324 


TOTAL 










-0.151 








-0.208 



Table 3: The three-body force contributions (in MeV) to the triton binding energy. The contribu- 
tions are arranged according to Eq. (20). Distinction between the isospin even and odd contribu- 
tions is also shown. The calculation was made perturbatively with a (three-channel) Malfliet-Tjon 
I-III model wave function. 
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FIGURE CAPTIONS 
Figure 1 — Diagram for the ir — n force 

Figure 2 — Diagrams for the tt — p force 

Figure 3 — Diagram for the p — p force 
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